Einstein— Yang— Mills Black Hole Interiors: 
Serious Problems But Simple Solution 
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Recently E. E. Donets, D. V. Gal'tsov, and the author reported the results of numerical and 
analytical in vestigation of the SU(2) Einstein- Yang-Mills (EYM) black hole interior solutions (gr- 
qc/9612067). It was shown that a generic interior solution develops a new type of an infinitely 
oscillating behavior with exponentially growing amplitude. Numerical data for three sequential 
oscillations were presented. The numerical integration technique was no t discussed. La ter P. Bre- 
itenlohner, G. Lavrelashvili, and D. Maison confirmed our main results (gr-qc/9703047). But they 



have made some misleading statements. In particular, they claimed, discussing the oscillations, that 
"as one performs the numerical integration one quickly runs into serious problems..." so that "it 
is practically impossible to follow more than one or two of them numerically" because "the nu- 
merical integration procedure breaks down" (pp. 3, 12). It is shown here that trivial logarithmic 
substitutions and integration along the integral curve solve these "serious problems" easily. 
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"Well", said Owl, "the customary 
procedure in such cases is as follows. " 

A. A. Milne ;) 



Recall the main notation from (T|. We assume the 
static spherically symmetric magnetic ansatz for the YM 
potential 

A = (W(r) - 1) (T v d6 - T e shi9d<p) , 

{Tipfi are spherical projections of the SU(2) generators), 
and the following parametrization of the metric: 

ds 2 = (A/r 2 )a 2 dt 2 - (r 2 /A)dr 2 - rW , 

where dil 2 = d9 2 +sin 2 9dip 2 , and A, a depend on r. The 
field equations include a coupled system for W, A: 



AU' + FW' = WV/r , 
(A/r)' + 2AU 2 — F , 



(1) 

(2) 



where U = W'/r, V = W 2 - 1, F = 1 - V 2 /r 2 , and a 
decoupled equation for a: 

(In a)' = 2rU 2 . 

These equations admit BH solutions in the domain r > r% 
for any radius t% of the event horizon. 

The system (JlJ, |J) was integrated numerically in the 
region < r < using an adaptive step size Runge- 
Kutta method starting at the left vicinity of the event 
horizon r/i with one free parameter Wh = W(rh) satis- 
fying inequalities \Wh\ < 1 and 1 — W% < r^, which are 
the necessary conditions for asymptotic flatness. 

One of the main results of our investigation is that a 
generic EYM BH interior solution develops a new type 



of an infinitely oscillating behavior with exponentially 
growing amplitude, during which the 'falls' of A (A < 0) 
may be approximated as 



A(r) 



A(r k ) 



r exp[U(r k ) 2 (r 2 - r 2 )} , r k > r , (3) 



where r k is the "starting point" of the k-th oscillation, 
{U{rk)r k ) 2 ^> 1, and all the regime can be described 
by a two-dimensional dynamical system (see [Q for the 
details). 

Our investigation of the EYM BH interior solutions 
|l]-|5| reveals, that one meets only two problems during 
numerical integration of (0, ||). 

Problem 1. Very small values of r as r — + and |A| in 
the neighborhoods of A's local maxima, and very large 
values of \U\ as r — > 0, and |A| in the neighborhoods of 
A's local minima. 

Problem 2. Intervals of very slow variation of r but 
extremely fast variation of A and maybe U . 

Obviously, the first problem can hardly be called a 
problem because it may be solved by trivial substitutions 
like ln(l/r) (or just lnr), lnln(l/r) , etc., for r; ln(— A), 
In | ln(— A)|, etc., for A an so on. 

The second one is not much more difficult. Let us show 
how it was solved for obtaining numerical data for three 
sequential oscillations, presented in JlJ (we shall also use 
figures 6, 7 from 01). 

Let us rewrite ([y, |^) in a more "comfortable" way, 
suggested in 0: 



DW" + FW' = —VW , 

x 



D' 



aW 



12 



2.r 



(4) 
(5) 



where x — {r/r^) 2 , D = 2A/i 



h> 



A/r 2 h , V = {W 2 



l)/2, F = 1 - aV 2 /x, and ' = d/dx. Here W' plays the 



1 



same role as U in ([l], ||). In a normal form, needed for 
the Runge-Kutta method, this system reads as 

W = P, 



P > = -Lv(W + aVP) - £ , 



(6) 



D' = F- 



aP 2 - 



1 

2x 



It is clear from (^|) that small (absolute) value of denom- 
inator xD, which takes place in the metric function local 
maximum neighborhood, may cause numerical problems. 
It may be unessential for some first maxima ('max 1' and 
'max 2' in Fig. [l]), but becomes considerable as oscilla- 
tions progress. Really, integrating system (0) for = 2, 
Wh = —0.342072, one finds out that x practically stops 
in the third local maximum neighborhood ('max 3' in 
Figs. |, |) with D » -10~ 16 , W w 10 18 (I mean a PC 
with 15 digits after the decimal point). This does not 
allow to pass the third maximum correctly. One of the 
possible solutions is to desingularize @ by introducing a 
parameter t as 

dx = xDdt. (7) 

But this method is in some sense local: it allows to pass 
the local maximum, but does not allow to reach the fol- 
lowing minimum of D. More general technique is based 
on integration along the integral curve. 

Recall the basic formula for the curve length (see, e.g., 
jjj or any other textbook on integral calculus). Let T be 
a smooth curve, defined by 



y 

where <fi, ip, x £ C^fa, b), and 
Then the length L of T is expressed as 

rb 

L 



t e [a, 6], 
X 12 > at 



VV 2 + V' 2 + X ,2 dt. 



This simple formula allows to solve our numerical prob- 
lems by introducing a parameter, representing the length 
of the integral curve, and to perform numerical integra- 
tion along the integral curve instead of the radial coordi- 
nate x (or r). This allows to pass the intervals of slowly 
varying x easily. 

Thus, one may introduce a parameter, say /, such 
that dl = (W 2 + P' 2 + D' 2 y/ 2 dx (or, for example, 
dl = (1 + P' 2 y/ 2 dx to integrate along P, etc.) and to 
rewrite system (j(^) in a relevant way. But we will intro- 
duce t = hix first in order to have the solution compo- 
nents of comparable values, and G = ln(— D) in order 
to pass not only the third local maximum, but also the 
third local minimum (see @ and Fig. ||). Thus, the new 
system may be written as 



(8) 



t 




- 1 
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dt 


e*P 
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~ ~dl 


NP - QVW 


. G _ 




_ | - ac l P 2 - N _ 




FIG. 1. The beginning of A-oscillations for ru = 2, 
W{r h ) = -0.342072 (n = 1 BH solution). 




FIG. 2. |A| for the second oscillation. 
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FIG. 3. ln(|A| _1 ) for the third oscillation. 

where Q = exp(-G), H = e* - aV 2 , N = HQ, ' = d/dl, 
and, for example, dl = (l + G 2 ) 1//2 <ii for integration along 
G. This system allows to pass the third local maximum 
and parts A and B of the third oscillation (Figs. S, 0). 
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FIG. 4. A versus I: the third local maximum. 



Note that part B in fact does not need to be passed along 
the integral curve. It is sufficient just to use t = lnx (or, 
e.g., t — ln(l/x)) as an independent variable. 

At last, one needs to pass part C (Fig. ||) in order 
to reach (and pass) the fourth local maximum of A. A 
vanishes fast at this interval; this forces W to increase 
strongly (\UA\ w \U(r™ ax )\ at this interval §). There- 
fore one needs to introduce Z = InW. System (||), 
rewritten in the corresponding way, solves this problem. 
(Evidently, this can be done before passing 'max 3', but 
(|J) works a bit better along part A.) 

It is also convenient to improve a system, prepared for 
integration along part C in a way similar to (]?]). Let 
us introduce another parameter, say p, such that dl = 
—Ddp. Then one obtains 



P 
i 

W 

z 

G 



dp 
~dl 



1 

E 

Eexp(t + Z) 

H - VW exp(-Z) 

\E -H -aE exp(i - 



■2Z) 



where E = exp(G). This system works fine both along 
part C and while passing local maxima (both 'max 3' and 
'max 4'). 

The next oscillation needs the next "order" of loga- 
rithmic substitutions (lnln(l/x), etc.), but the numerical 
integration technique remains the same. 

Recall, at last, that the investigation, presented in [Q, 
was restricted to the class of asymptotically flat solutions. 
If one also studies solutions, which are singular in the 
exterior region, then a situation, shown in Fig. ||is rather 
typical. This makes the claims cited in the abstract even 
more surprising. 

Thus, if we do not discuss restrictions, which are in- 
herent for every numerical method, there are hardly any 
"serious problems" , which can "break down" the numer- 
ical integration of the SU(2) EYM equations inside the 
event horizon. All the difficulties, connected with the 
numerical investigation, are typical and can be solved by 
integration along the integral curve. 
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FIG. 5. This curve represents a solution, singular in the 
exterior region. One can see 4 first oscillations of A and the 
beginning of the 5th one, which can be easily integrated using 
the technique, described above. 

I would like to thank Rostislav Zhukov for his little 286 
PC. The research was supported in part by the RFBR 
grant 96-02-18899. 
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